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Abstract 



We make for the first time a large-scale Monte-Carlo simulation of a fer- 
romagnetic Heisenberg model with dipolar interactions on a two dimensional 
square lattice with open boundaries using an efficient new technique. We 
find that a phase transition occurs in the model and the ordered phase is 

characterized by a circular arrangement of the spins. 
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Spin structures in systems with competing interactions have been of great interests in 
the last decades. In the homogeneous systems, competition between the nearest neighbor 
and the second nearest neighbor interactions brings a hehcal spin ordering in the Heisenberg 
model and a devil's staircase in the Ising model 0. In the random systems, spin glass 
phase occurs due to the competition between ferromagnetic and antiferromagnetic exchange 
interactions 1^,^ and the competition in the anisotropy between two ions in the antifer- 
romagnetic mixture causes a new magnetic ordering called as an oblique antiferromagnet 
In these cases, magnitudes of the competing interactions are of the same order and 
the spin structures are modulated microscopically. Another kind of competition occurs in 
the ferromagnetic system, namely, the exchange interactions compete with the dipolar in- 
teractions which are much smaller than the formers. This competition has been studied in 
technological interest rather than the physical one, because it realizes magnetic domains as 
in the real materials. Magnetic structures of ferromagnetic thin films have been attracted 
in recent years, especially in the field of the magnetic storage technologies p|-p!0|]. How- 



ever, the models studied so far are phenomenological ones in which systems are divided into 
small cells with a certain magnetization and anisotropy energy and the distribution of the 
magnetizations is determined using the Landau-Lifshitz-Gilbert equation [^T]-|T^. 

Although the models reproduce the domain structures of those films fairly well, there 
still remains primary questions whether or not those structures could really be reproduced 
in a microscopic model and, if it could, whether or not a clear-cut phase transition could be 
seen? To answer the questions, we are necessary to make a large-scale Monte Carlo (MC) 
simulation of a Heisenberg model with dipolar interactions, because the domains can be seen 
only in large systems. The bottleneck of the simulation is, however, to calculate the dipolar 
fields of the individual spins which needs the CPU time of order A^^ in every MC step, where 
A^ is the number of the total spins. Because of this reason, no large-scale MC simulation 
study has been done up to now on the spin structure as well as the phase transition of this 
model |TB|,|TB . 



In this letter, we report for the first time results of a large-scale MC simulation of a 



ferromagnetic Heisenberg model on the two dimensional square lattice with open boundaries. 
We develop an efficient new MC technique which reduces the CPU time of order from N"^ to 
N log2 N . We find that a clear-cut phase transition really occurs at a finite temperature and 
a domain-like spin structure appears below that temperature. The structure is shown to be 
characterized not by closure ferromagnetic domains but by a circular arrangement of the 
spins. We call, hence, the ordered phase as a circular phase. Values of the critical exponents 
are also estimated to examine the nature of the phase transition. 

We start with the classical Heisenberg model described by the Hamiltonian: 

<i,j> ijtj ^ij ^ij 

where = 1. The former term describes the exchange energy, where J(> 0) denotes the 
exchange constant and runs over the nearest neighbor pairs. The latter term describes 
the dipolar energy, where D — [q^bSY /a^ but being considered here as a parameter and 
Tij is the distance between ith and jth spins which is measured in the units of the lattice 
constant a. Hereafter we measure J and D in the units of /c^ = 1 with fc^ being the 
Boltzmann constant. The lattice is the square lattice oi L x L with open boundaries. We 
choose the x- and y-axes being parallel to the lattice plane and the 2;-axis perpendicular to 
them. 

The MC simulation is made as follows. Since the dipolar interactions are of long-range, 
huge CPU time is necessary to calculate the dipolar field H^^ for every spin 5,. Then we 
distinguish the dipolar field H^^ into two parts and i?^2> i-^-> contributions 
of the spins up to the second nearest neighbors and Hj^2 those of the other spins. We 
calculate H^^i well as the exchange field for every update trial. On the other hand, H^o2 
is updated for every m MC steps(here we choose m = 5). Therefore, we can calculate the 
set of the dipole fields {H^d2} simultaneously by using FFT algorithm which needs the CPU 
time of order A^" log2 N. Of course, errors arise which will become larger as the temperature 
is decreased. However, we think that the errors are not so serious in the thermal equilibrium 
because of the following reasons. Since a large number of the spins may not change their 
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directions in each MC step, the change of the dipolar field H^]j2 due to the changes of far spins 
will not be accumulated considerably. Then the major part of the errors will come from the 
ignorance of the updates of near spins, especially those on the third nearest neighbor lattice 
sites. Therefore, the errors will become considerable below T~4xD/2'^ = D/2. However, 
at such low temperatures, the spins will change their directions very little in each MC step. 
Hence we believe our method does not change results significantly at all temperatures ||17|| . 
We carry out the simulation of the model with J = 1 and = 0.1 using the heat bath 



method [|18|. The linear sizes of the lattice are L = 16 — 256. We make the MC simulation 
for several times starting from different initial spin configurations. The MC steps are about 
100,000 for the largest lattice. Coohng the temperature, we calculate the specific heat C, 
the absolute value of the uniform magnetization and the susceptibility (/^ = 2^,1/, z): 

M„ = < \M\ >, (3) 
X" . SiM^, (4) 

where M = {l/N)J2i Si, and < ... > means the thermal average. We also observe the spin 
structure. 

In Fig. 1, we show typical three spin structures A, B and C observed at low tempera- 
tures. The structure A is single-domained and the structure C like domains of closure which 
frequently observed in thin films 0. The structure A is seen for smaller lattices of L < 16 
and the structure C for larger lattices of L > 32. The structure B is an incomplete one 
which occurs for intermediate size of L ~ 24. This relation between the spin structure and 
the lattice size is not so clear-cut. For L = 24, the structure B dominates over the structures 
A and C, but for L = 32 the structure C is observed most frequently. 

We show the specific heat C and the z-component of the susceptibility in Figs. 2 and 3, 
respectively. The specific heat exhibits a size dependent sharp peak at T ~ 0.85. Hereafter, 
we denote the peak temperature as Tc. The susceptibility also exhibits a size dependent 
maximum at almost the same temperature of T^. These results suggest the occurrence of 



some long-range order characterized by the xy-components of the spins. If the long-range 
order really occurs, it would be one with the spin structure C shown in Fig.l, because we 
could not observe any evidence of another phase transition below Tc. Then we calculate, in 
addition to M^, the absolute value of the circular component of the magnetization M^: 



where and rc are the position vectors of the ith spin and the center of the lattice, re- 
spectively. In Fig. 4, we show temperature dependences of M„ for different L. When L is 
small, Mu increases as the temperature is decreased and seems to reach 1 at T = 0. How- 
ever, as L is increased, M„ only exhibits a hump around T^. In Fig. 5, we show temperature 
dependences of for different L. We see an opposite size dependence, i.e., for smaller 
L exhibits a hump around Tc, while for larger L exhibits a rapid increase below Tc. 
Moreover, the size dependence of is reversed as the temperature is decreased below Tc. 
That is, for a larger L has a larger value below T^. These results clearly reveal that the 
long-range order really occurs and it is described not by the uniform magnetization but by 
the circular component of the magnetization. 

We should mention that the ordered phase obtained above is not the one composed of 
four ferromagnetic domains with their magnetization axes parallel to the edges of the lattice, 
because if so, one would get ~ sinh~^(l) ~ 0.88 at T = 0, whereas the result in Fig. 5 
suggests ~ 1 at T = 0. Moreover, we could not see any domain-like structure near below 
Tc. Therefore it is natural to call the phase as a circular phase rather than the ferromagnetic 
phase with domains. Next we examine the phase transition of the model using the finite- 
size-scaling analysis. If the second order phase transition occurs at T = Tc, for the 
lattice with L will be scaled as 



where t = {T — Tc)/Tc, (3 and v are exponents of the order parameter and the correlation 
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Ti - r. 



(5) 




(6) 



length, and F is some scaling function [^. We could fit the data for L > 64 fairly well, 
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which is shown in Fig. 6. From the results, we estimate values of the exponents as /3 ~ 0.23 
and V ~ 1.20 together with ~ 0.88. We have also estimated the values of (5 and a by 
plotting vs T and C vs T in log-log forms, respectively, and obtained similar values of 
/3 ~ 0.2, a ~ —0.4 and Tc ~ 0.87. Note that this value of a is compatible with the scaling 
relation oi a = 2 — dv with d = 2. These values are not similar neither to those of the Ising 
model of /3 = 1/8 and u = 1 nor to those of an anisotropic xy-like model without dipolar 
interactions ||20|. We conclude, hence, that a usual second order phase transition really 
occurs in this model and it belongs to a new universality class. This is another evidence 
that the circular phase is different from the ferromagnetic phase with domains. We note 
here that the present result is not incompatible with experimental observations. If an axial 
anisotropy along the edges of the lattice exists, the spins tend to align the directions parallel 
to the edges. Thus the circular structure is deformed into the one with four ferromagnetic 
domains of closure which are observed in thin films [§. Finally we should emphasize that 
the results presented here have been obtained by using the new MC technique within the 
CPU time much less than that of the conventional MC method. 

The authors would like to give their thank to Dr. K. Tan for leading them to this problem, 
and to Professor T. Shirakura and Dr. T. Nakamura for their valuable discussions. This 
work was financed by a Grant-in-Aid for Scientific Research from the Ministry of Education, 
Science and Culture. 
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FIGURES 

FIG. 1. Snapshots of the spin structures for different sizes of the lattice at T = 0.1. 
FIG. 2. Temperature dependences of the specific heat C for different sizes of the lattice. 

FIG. 3. Temperature dependences of the susceptibility for different sizes of the lattice. 

FIG. 4. Temperature dependences of the absolute magnetization for different sizes of the 
lattice. 

FIG. 5. Temperature dependences of the circular component of the magnetization for 
different sizes of the lattice. 

FIG. 6. A finite-size-scaling plot of M^. 
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Fig.3 
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